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Abstract. Earthworm distributions are poorly known and individual species climatic tolerances, even less so. This paper sets 
out to use three species with a mainly Anglo-French distribution to test out whether using Maximum Entropy Modelling 
(Maxent) could be useful when studying earthworm distributions. It also gives an indication of how the likely climatic 
changes over a 50 year period will affect them. Overall the software seems to give useful information of where across Europe 
a particular species will thrive, even if not currently recorded there. It gives a real insight into how particular species might be 
better able to survive longer drier periods than others and which are on the edge of their climatic range already. Maxent 
modelling was clearly successful in demonstrating that the distributions of the ecologically different earthworm species are 
affected by a combination of different environmental variables. In the case of the epigeic Satchellius mammalis they are the 
annual temperature range, the precipitation of the driest month and the mean annual precipitation, for the epi-endogeic 
Lumbricus festivus they are the precipitation of the driest month, the precipitation of the wettest month and the annual 
temperature range. For the anecic Lumbricus friendi the most important environmental variables proved to be the annual 


temperature range, the mean diurnal temperature range and the precipitation seasonality. 
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INTRODUCTION 


[ is well known that different earthworm spe- 
cies have very different environmental toler- 
ances (Lee 1985). For example some earthworm 
species have cocoons that are more frost resistant, 
others more drought resistant (Holmstrup 1994). 
To predict the possible impact of climate change 
software such as Maxent (Phillips et al. 2006) 
developed by the American Museum of Natural 
History, can be used. This software allows the 
input of detailed maps of environmental variables 
as well as the location data of specimens. Using 
Maxent, and suitable environmental variables, it is 
also possible to model the effects of future cli- 
mactic conditions on the suitability of different 
environmental locations for the different earth- 
worm species. Consequently, potential future pre- 
dictions with regards to habitat change or climate 
change could be considered. 


Maxent was successfully used to predict the 
distribution of an endogeic earthworm Hormo- 
gaster elisae Alvarez, 1971 in central Spain using 
climatic and soil variables (Machan et al. 2015). 
Latif et al. (2017) used Maxent to model the dis- 
tribution of the epigeic sibling species pair Eise- 
nia fetida (Savigny, 1826) / Eisenia andrei Bou- 
ché, 1972 in Iran and showed that the most 
important environmental variables in determining 
the natural distribution of E. fetida/andrei were 
annual mean temperature and precipitation in the 
driest months followed by the mean diurnal range 
of temperature and precipitation in the wettest 
months. Also, using Maxent modelling Geraskina 
& Shevchenko (2019) successfully demonstrated 
that the main climatic factor influencing the 
distribution of the two epigeic species Dendro- 
baena octaedra (Savigny, 1826) and D. attemsi 
(Michaelsen, 1902) in the northwestern Caucasus 
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is the precipitation in the driest month which has 
an effect on the desiccation of the litter layer and 
so on the survival of the epigeic worms in the 
driest summer months. 


With ever increasing understanding of the im- 
portance of earthworms as ecosystem engineers 
(Lavelle et al. 2016), understanding the likely 
effect of climate change on distribution of indi- 
vidual species is crucial. However, testing of 
environment modelling like Maxent on conti- 
nental scale data is still very rare, mainly due to 
the lack of continent-wide datasets. 


Earthworm distribution globally is quite poorly 
known (Blakemore 2010). The UK, despite 
having a relatively small fauna with full iden- 
tification keys present since the late 1940’s (Sims 
& Gerard 1985), mirrors this poor distributional 
data (Carpenter et al. 2011). With the estab- 
lishment of the earthworm society of Britain and 
in particular its launch of NERS (National Earth- 
worm Recording Scheme) this situation is gra- 
dually starting to change now with over 12,000 
good quality data records which together with 
other available datasets are analysed here. 


We hope this preliminary analysis can be the 
beginning of future detailed continental-scale 
work making informed judgments on how cli- 
matic changes and habitat destruction might be 
affecting these important soil ecosystem engineer 
taxa. 


MATERIALS AND METHODS 


In this study, Maximum Entropy Modelling 
(Maxent?) (Phillips et al. 2006) was used to 
investigate the potential distributions of different 
earthworm species. Using the specimen location 
data, as well as certain bioclimatic data that we 
provided, Maxent produced a map detailing the 
suitability of different environmental locations for 
each species. 


Specimens. The data for the species used were 


collected in three different ways. The first was 
through a compilation of Museum collections data 
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and private research data. Specifically these were; 
The Natural History Museum London with 
collections data from the late 1800’s to the present 
day (400 relevant records), The Hungarian Na- 
tural History Museum Budapest with collections 
dating from the middle 1900’s to the present day 
(approx. 80 relevant records), The Smithsonian 
Institution in Washington DC with collections 
data from the late 1800’s to the present day 
(approx. 25 relevant records) and 8 records of 
Satchellius mammalis from Sweden and Norway 
from private research data held by Christer Erseus 
in Sweden collected between 2008 and 2012. 


For each piece of collections data which con- 
tained a vague location, an extensive search was 
carried out to find that location and obtain its 
latitude and longitude. This was not always pos- 
sible. This may have been due to place name 
changes, or places with the same name and region 
such that the location given was not specific 
enough to differentiate, or errors in recording the 
location given. In such cases, that data record was 
disregarded. 


Alongside collections data, data records from 
the Earthworm Society of Britain were also used 
to provide detailed information about earthworm 
populations in Great Britain. This data includes 
some historical records but most records were 
collected in the past 10 years (approx. 10,000 
records). The location of sample, as well as lati- 
tude and longitude are recorded. 


The third way data was collected using the 
book of Bouché (1972) where he recorded large 
amounts of location data for different earthworm 
species across France. Although this data is now 
50 years old the climatic conditions are mapped to 
the time recorded at so the data is still very valid. 
The data for the each species used in this study 
were compiled and latitudes and longitudes for 
the locations were searched via Google Maps. 


Three species of earthworms were selected for 
detailed mapping across Europe; Satchellius 
mammalis (Savigny, 1826), Lumbricus friendi 
(Cognetti, 1904) and Lumbricus festivus (Savigny, 


Sherlock, Coates & Csuzdi: Modelling of climatic tolerances of three earthworm species using Maxent 


1826). These species were selected because their 
restricted ranges are centering on the Anglo- 
French region from which the greatest proportion 
of our records are, and because these species had 
limited ranges therefore it could be deduced that 
they are particularly sensitive to changes in the 
climatic variables. 


For investigating the effect of the land use on 
the range of widely distributed peregrine earth- 
worms two of the UK’s most common earthworm 
species were selected; the endogeic Aporrectodea 
caliginosa (Savigny, 1826) and the epi-endogeic 
Lumbricus rubellus Hoffmeister, 1843. 


Due to the possibility of sampling bias in re- 
gards to the Earthworm Society of Britain, a map 
containing the location of every record of the so- 
ciety was also created. 


The model. The model was created using 
eleven different bioclimactic variables from 
WorldClim 1.4, at a size of 30 arc-seconds. This 
was to allow high detail on smaller locations, such 
as the British Isles. A selection of important 
variables relating to both temperature and preci- 
pitation levels was chosen. The bioclimatic vari- 
ables chosen were: 


1. Annual Mean Temperature 
This variable shows the mean temperature of a 
location for a single year. 

2. Mean Diurnal Range 
This variable demonstrates the mean range in 
temperatures on a single day. 

3. Maximum Temperature of the Warmest Month 
This variable demonstrates the highest tempe- 
rature that occurs in a location, in the month 
which is on average warmest in that location. 

4. Minimum Temperature of the Coldest Month 
This variable demonstrates the lowest tem- 
perature that occurs in a location, in the month 
which is on average coldest in that location. 

5. Temperature Range 
This variable shows the range in temperature 
between the average temperature of the 
warmest month and the average temperature 
of the coldest month. 

6. Mean Temperature of the Wettest Quarter. 


This variable demonstrates the mean tem- 
perature during the wettest three months of the 
year in that location. 

7. Mean Temperature of the Driest Quarters 
This variable shows the mean temperature 
during the driest three months of the year in 
that location. 

8. Annual Precipitation 
This variable shows the amount of preci- 
pitation a location gets over an entire year. 

9. Precipitation of the Wettest Month 
This variable shows the amount of precipi- 
tation a location gets during the month with 
the most precipitation. 

10. Precipitation of the Driest Month 
This variable shows the amount of precipita- 
tion a location gets during the month with the 
most precipitation. 

11. Precipitation Seasonality 
Precipitation seasonality measures the varia- 
tion of precipitation totals between each 
month over the year. 


The temperature variables were chosen to re- 
flect the effect of temperature fluctuations on the 
distribution of different earthworm species on dif- 
ferent time scales. For example, Mean Diurnal 
Range was chosen to reflect how daily tem- 
perature fluctuations may affect earthworm dis- 
tribution, whilst Temperature Range and Annual 
Mean Temperature were chosen to reflect how 
temperature fluctuations over an entire year may 
affect earthworm distribution. Similarly Maxi- 
mum Temperature of the Warmest Month and 
Minimum Temperature of the Coldest Month 
were both chosen to investigate the possibility 
that some earthworms may prefer mild climates, 
and the effects of high and low temperatures on 
earthworm activity and the durability of their 
cocoons. For example, in the case of Minimum 
Temperature of the Coldest Month, this may be 
due to soil freezing. 


The precipitation variables were chosen to 
reflect how the availability of moisture in the soil 
fluctuates through the year. Annual Precipitation 
was chosen to reflect how much water is available 
in a location during a year. The Precipitation of 
Wettest and Driest months was chosen to reflect 
how a sustained period of large amounts of water, 
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or lack of water, in the soil may affect earthworm 
populations, especially those not known to aes- 
tivate in the driest months. Precipitation Season- 
ality was chosen to reflect how large or small va- 
riations in precipitation over the year may affect 
earthworm populations. 


Future projections of each of these bioclimatic 
variables were also used to produce the future 
projections of the suitability of different environ- 
mental locations for each species. Specifically we 
used the RCP 6.0 (Representative Concentration 
Pathway) pathway projected to 2050. The reason 
we chose RCP 6.0 is because it was the higher of 
the two middle RCP’s, providing a possible worst 
case scenario, without using RCP 8.5 which may 
be overestimating future supply of fossil fuels 
(Rutledge 2011, Wang et al. 2017). We also did 
not attempt to produce maps for climate pro- 
jections after 2050, due to that providing too 
much uncertainty. 


The importance of different predictor variables 
for each species analyzed was determined accord 


ing to Analysis of Variable Contributions (AVC). 
The variables that make a significant contribution 
to the model are those which have high values of 
permutation importance (PI) (Phillips et al. 2006). 


In the case of Great Britain, land use was also 
investigated to ascertain whether it also has a 
large effect on earthworm distributions. Land use 
maps were taken from the Centre for Ecology and 
Hydrology (CEH). Land use for Ireland and 
Northern Ireland were not included. 


RESULTS 


1. Satchellius mammalis shows mainly Anglo 
French distribution which is primarily concen- 
trated around the Western North coast of France, 
the Netherlands and Belgium. According to Max- 
ent, all regions of the UK, except exposed wester- 
ly areas of Northern Ireland, Westerly highland 
areas of Scotland and an area of eastern England, 
are favourable. Hotspots seem to be centred a- 
round South coastal regions and the Welsh Eng- 
lish border (fig. 1). 


Figure 1. Satchellius mammalis. Present European distribution map with climatic variables. 
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Figure 2. Satchellius mammalis. Predicted future European distribution map with climatic variables. 


Table 1. Importance of different climatic variables in determining the distribution of S. mammalis according to the AVC data. 


Variable Percent contribution | Permutation importance 
Temp Annual Range 39.7 43.7 

Precip Driest Month 22 27.4 

Mean Precipitation 1.9 9.7 

Mean Temp 0.7 6.6 

Mean Diurnal Range 1.5 5.1 

Mean Temp Wettest Quarter 3.6 3 
Precipitation Seasonality 4 2 
Min Temp Coldest Month 25.2 1.3 
Mean Temp Driest Quarter 0.3 1.2 
Precip Wettest Month 1.1 0 

Max Temp Warmest Month 0 0 


The main variables determining the distri- 
butions according to AVC (Table 1) are: tem- 
perature annual range and the precipitation in the 
driest month. This would suggest this species is 
less drought resistant. It fits well to the maps 
produced as the area of East Anglia is one of the 


driest in the country and this species seems to fare 
less well in those drier more easterly regions 


Modelling the future distribution of S. mam- 


malis shows a predicted widening of range in the 
more northerly regions it inhabits but the souther- 
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ly regions become less favourable (fig. 2). This is 
especially true in the UK with a predicted larger 
scale increase in optimal conditions like higher 
annual temperature in the northern regions and 
more precipitation in the summer. 


If we incorporate in the model the UK land use 
map it is clear that the extent of the most favour- 
able area for this species notably shrinks but its 
geographic location in the UK still remains 
largely the same (fig. 3). 


Using the climatic predictions set in combi- 
nation with the land use data a possible south 
western shift of the species’ favourable range 
would appear in comparison with the modelling 
without land use data (fig. 4) 


2. Lumbricus friendi distribution (fig. 5) shows 
that the UK is only on the very fringes of this 
species climatic range currently. South eastern 


France, especially around the mountainous areas 
such as the Pyrenees and Massif Central, are 
favourable. A favourability towards mountainous 
regions is also supported by a clustering around 
the Alps in the East. However, they are not found 
in the highest altitudes of these ranges, just the 
surrounding areas. The most important predictors 
of the present distribution according AVC (Table 
2) are; the temperature annual range, mean diurnal 
temperature range and precipitation seasonality. 
Future predictions with the warmer climatic 
conditions do not appear to be favourable for L. 
friendi (fig. 6) as the most favourable areas dis- 
appear almost entirely in the predictions, however 
the extent of the suboptimal areas especially in 
northern France and in the Ardennes, Belgium 
seems to be increasing. Taking into account the 
semi-peregrine nature of this species (Csuzdi & 
Szlavecz, 2004) this might indicate a possible 
North-Eastward shift in its distribution. 


Figure 3. Satchellius mammalis. Present UK distribution 
map with the addition of land use data. 
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Figure 4. Satchellius mammalis. Predicted future UK 
distribution map with land use data. 
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Figure 6. Lumbricus friendi. Predicted future European distribution map with climatic variables. 
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Table 2. Importance of different climatic variables in determining the distribution of L. friendi according to the AVC data. 


Variable Percent contribution | Permutation importance 
Temp Annual Range 23.1 31.9 

Mean Diurnal Range 22.3 24.2 
Precipitation Seasonality 27.3 15.8 
Mean Temp 3.9 7.7 

Precip Driest Month 6.9 6.5 

Precip Wettest Month 4.3 5 

Mean Temp Wettest Quarter 4.8 3.6 
Min Temp Coldest Month 0.4 2.3 
Mean Temp Driest Quarter 5.8 1.7 
Mean Precipitation 0.4 0.9 

Max Temp Warmest Month 1 0.4 


3. Lumbricus festivus shows a typical Atlantic 
distribution type (Csuzdi & Zicsi 2003) centring 
in Southern and Midland regions of the UK and 
the more northerly France. The major areas of 
unsuitability in the UK are the areas of higher 
altitude such as the Lake District, Pennines, 
Western Scotland and central Wales. A large area 
of eastern England is also less favourable (fig. 7). 


The predicted favourable locations for this 
species remain mostly the same with the future 
predictions with indication of a slight Northern 
shift of its range in France and also in UK. This is 
especially prominent in the UK, where most of 
the country becomes incredibly favourable, in- 
cluding the previously unsuitable areas of eastern 
England (fig. 8). The AVC (Table 3) indicates 
that, like S. mammalis, this species is not very 
drought resistant. The driest areas in the UK, in 
the east, are areas the species is not thought to be 
found. The second largest factor is precipitation in 
the wettest month then the third is temperature 
annual range. So this species is very sensitive to 
the amount of rainfall the region has and therefore 
is unlikely to survive prolonged dry periods. 


Incorporating the land-use data into our model 
shows a negative effect on the predictions, 
especially in a large area of Eastern England 
which becomes more prominently unsuitable (fig. 
9). This restricted range largely overlaps with the 
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present distribution of the species and explains its 
lack from Eastern England instead of the illus- 
trated favourable climatic conditions. 


In the future prediction, as with fig. 9, larger 
areas of the UK become more favourable for the 
species (fig. 10). However, the increased favoura- 
bility is more concentrated on the southern and 
north westerly coasts of England and the Eastern 
regions still remain unfavourable. This shows that 
land use data does have a significant negative 
effect on the predictions for this species, com- 
pared to fig. 8. 


To demonstrate if land use has an effect on 
widespread species we have modelled the dis- 
tribution of Ap. caliginosa (figs. 11-12) and L. 
rubellus (figs. 13-14) two of the most common 
species throughout Europe. Apart from exposed 
Westerly regions these species are common in 
most areas. According to the model predictions 
incorporation of land-use data makes almost no 
difference to the range of these species. 


DISCUSSION 


The AVC tables produced by Maxent analysis 
proved to be very useful in determining the spe- 
cific climatic conditions that each species is most 
susceptible to. For two of the three species 
analysed (S. mammalis and L. festivus) precipi 


Sherlock, Coates & Csuzdi: Modelling of climatic tolerances of three earthworm species using Maxent 


Figure 8. Lumbricus festivus. Predicted future European distribution map with climatic variables. 
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Table 3. Importance of different climatic variables in determining the distribution of L. festivus according to the AVC data. 


Variable Percent contribution | Permutation importance 
Precip Driest Month 14.6 33.2 

Precip Wettest Month 3.6 18.8 

Temp Annual Range 41.9 17 
Precipitation Seasonality 6.9 9.3 
Min Temp Coldest Month 16.8 8.5 
Mean Temp Driest Quarter 9.2 6.8 
Mean Temp 1 2.8 

Mean Temp Wettest Quarter 3.5 2 
Mean Diurnal Range 0.9 1.4 

Mean Precipitation 0.9 0.2 

Max Temp Warmest Month 0.7 0 


Figure 9. Lumbricus festivus. Present UK distribution 
map with land use included. 


tations in the driest months were among the most 
significant predicting factors, indicating that they 
are not drought tolerant. However, the third 
species (L. friendi) was more susceptible to tem- 
perature ranges. 


Bouché (1972) in his extensive volume on the 
lumbricid worms of France presented an excellent 
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Figure 10. Lumbricus festivus. Predicted future UK 
distribution map with land use included. 


insight into the climatic tolerances of the 
individual species including Satchellius mam- 
malis. It is a species with a healthy distribution in 
the UK and France, with the western and northern 
French coast and the southern British coast being 
particular hotspots. Bouché (1972) stated the 
species was likely to be restricted to the Northern, 
Western and central areas of France, not being 
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found below a line running from Strasbourg to 
Perpignan. He also stated however, that if intro- 
duced there, they would also thrive in a small area 
of Italy. When modelling the data available, 
Maxent predicted exactly that (fig. 1), so illu- 
minating the incredible insight by Bouché but 
importantly also giving weight to the reliability of 
Maxent to correctly predict species tolerances 
even from a relatively small sample size. 


One interesting area to explore in the UK for 
this species is a large area not considered habi- 
table around Peterborough in the South, Leicester 
in the west, Grantham in the north and Boston in 
the east (fig. 3). This needs a lot more investi- 
gative work as to exactly what the factors are 
affecting this. We assume, especially with an epi- 
geic species like S. mammalis,, which is not able 
to aestivate in the summer, the most likely 
explanation would be due to this area having the 
driest conditions in the UK. When the potential 
future predictions are considered it seems that this 
area will become more habitable, which again 
could be due to a higher predicted rainfall due to 
climate change (fig. 4). 


We considered under sampling could be skew- 
ing the presumed tolerances for this species. It is 
possible that the Earthworm Society of Britain 
does not sample as much in the east of England. 
To investigate this we mapped the location of 
every sample collected by the Earthworm Society 
of Britain (fig. 15). As can be seen on the map, 
there does appear to be a degree of under sam- 
pling in the east of England. However, we do 
have samples there so this alone does not explain 
Maxent showing a lack of suitable environmental 
conditions in this area. If the east of England 
contained suitable environmental conditions for 
the species, it is likely that Maxent would have 
predicted this based on the extensive sampling 
from the rest of Great Britain. 


One of the rare lumbricid earthworms in Eu- 
rope is Lumbricus friendi. There are just four 
records for the UK known of this species but we 
have numerous records from France, where the 


main extent of its worldwide distribution can be 
found. The records however, appear only to be 
found in very particular habitat types, and it 
would seem associated with certain (not high) 
altitudes, and the cooler moister climes which 
characterise them (fig. 5). This is demonstrated by 
a clustering of the species around the bases of the 
Pyrenees, the Massif Central, and the Alps, 
however, not being found higher in the mountain 
ranges themselves. Altitude though cannot be the 
sole contributor to this distribution, with speci- 
mens also being found in lower lying areas further 
away from these mountain ranges. This species 
seems to be particularly sensitive to climate with a 
very limited range of tolerable conditions, as 
shown by the future predictions largely reducing 
its suitable habitats (fig. 6). The climatic variables 
which are most important to this particular species 
seem to be precipitation seasonality and annual 
temperature range. So it would appear they need a 
constant rainfall with limited variation between 
months and a low but not freezing temperature but 
this need a lot more investigation. 


The Maxent map for L. friendi, a species con- 
sidered recently for a red list in the UK, shows it 
is absolutely on the boundary of its climatic tol- 
erances in the UK and is not a species that is 
suffering due to any particular habitat loss or 
change but a species on the absolute fringes of its 
ecological tolerance (and not likely to be much 
aided from being added to a red list) and never 
likely to thrive. With a species with such narrow 
tolerances the future could be disastrous, however 
the 50 years predictions again seem to be sug- 
gesting a more rosy outlook for the species in the 
UK albeit not anywhere else in its range (fig. 6). 


It is interesting that modelling identified an- 
other highly suitable region for L. friendi along 
the Adriato-Mediterranean region. Just recently, 
Stojanović et al. (2014) reported this species for 
the first time from Serbia. Although this record 
may be due to introduction like the North Ameri- 
can ones (Csuzdi and Szlavecz, 2004) it also de- 
monstrates the predicting power of Maxent ana- 
lysis with its highlighting of favourable regions. 
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Figure 11. Aporrectodea caliginosa. Present UK Figure 12. Aporrectodea caliginosa. Present UK 
distribution map without landuse information. distribution map with landuse information. 


Figure 13. Lumbricus rubellus. Present UK distribution Figure 14. Lumbricus rubellus. Present UK distribution 
map without landuse information. map with landuse information. 
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Figure 15. Map of the British Isles showing the Earthworm 
society of Britain’s sampling localities. 


Another species of interest is Lumbricus fes- 
tivus (fig. 7). This species is not as rare as L. 
friendi however, is still considered rare in the UK. 
Within France it appears to have a northern dis- 
tribution which is also predicted will extend out 
through Belgium and into the Netherlands. This 
would need to be tested via records from these 
countries. Similarly in the UK this species has a 
southerly distribution. It consequently appears 
that the areas around the English Channel and 
southerly North Sea seem to provide the most 
favourable conditions for this species. However, 
Maxent does also predict a very favourable chan- 
nel through France running from Le Havre and 
the surrounding area right through to Lyon and 
beyond. This again would need to be explored 
further as this does not quite match the records 
produced by Bouché (1972). 


But how much could land use play a role in the 
distribution of earthworm species? Another useful 
feature of Maxent is that other maps can be over- 
layed as well, so we used in this paper CEH land 
use map. Looking at very common Europe wide 


species with broad climatic tolerances such as 
Aporrectodea caliginosa (figs. 11-12) and Lum- 
bricus rubellus (figs. 13-14) we set to test if land- 
use were to affect their range. The results 
appeared to show there is very little difference 
between maps solely based on climate and when 
land use was added for these species. These 
worms can survive in most areas whether it be 
urban, agricultural or wild. However, when com- 
paring this to a map of Satchellius mammalis 
(figs. 1, 3), one of the more ’sensitive’ species, we 
can see that in this case land use has had a clear 
effect. Areas in the east in particular suggest a 
dearth of records even beyond that predicted for 
the dry climatic conditions alone, which our 
current data set supports. Looking at the land-use 
maps closely this seems to tie in with this species 
struggling in areas being converted to arable 
lands. However, agriculture based more on pas- 
tures for livestock do not have such a significant 
effect, as shown in the more westerly regions. So 
although this feature needs to be investigated 
further, in this limited study land use does seem to 
be a useful factor to consider for species with 
limited ranges. As we do not currently have de- 
tailed land use maps for Europe it might also be 
useful to investigate the effect of land use on the 
favourability of regions in France for certain 
species. This may be especially enlightening in 
the case Lumbricus friendi, which we know to be 
sensitive to climactic conditions, and may also be 
sensitive to different land uses. 


CONCLUSION 


This is an initial study to test Maximum En- 
tropy Modelling on a limited number of earth- 
worm records and species. Species distribution 
models of course always have limitations (Jar- 
nivich & Young 2015, Carneiro et al. 2016) 
however, the initial results do seem to indicate 
that Maxent is a useful tool to use in ascertaining 
general climatic tolerances and trends in the 
distribution of individual earthworm species; and 
that it would be worthwhile investing time and 
energy in a broader study with a greater number 
of records with a much broader geographical 
distribution (Rutgers et al. 2016) 
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This was especially evident with the case of L. 
friendi in the west Balkans and in France. It seems 
in particular to be useful with those species unable 
to cope with large fluctuations of temperature, or 
moisture. It therefore could also, when future pre- 
dictions can be made of precipitation and tem- 
perature fluctuations enable us to broadly predict 
the likely vitality and range of these species. 


Our results clearly emphasises that a lot more 
work needs to be carried out in collating all 
known verified records of European earthworm 
species as well as systematically collecting large 
quantities of species level data. From this study 
however, we do feel, once generated, this infor- 
mation alongside the use of environmental model- 
ing like Maxent, should start to help us have a 
much healthier picture of how our earthworms are 
faring and how they are likely to fare in the 
future. 
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